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Two-dimensional dilute magnetic semiconductors can provide fundamental insights in the very 
nature of magnetic orders and their manipulation through electron and hole doping. Despite the 
fundamental physics, due to the large charge density control capability in these materials, they can 
be extremely important in spintronics applications such as spin valve and spin-based transistors. In 
this article, we studied a two-dimensional dilute magnetic semiconductors consisting of phosphorene 
monolayer doped with cobalt atoms in substitutional and interstitial defects. We show that these 
defects can be stabilized and are electrically active. Furthermore, by including holes or electrons 
by a potential gate, the exchange interaction and magnetic order can be engineered, and may even 
induce a ferromagnetic-to-antiferromagnetic phase transition in p-doped phosphorene. 


INTRODUCTION 

Over the last decades the field of dilute magnetic semi¬ 
conductors (DMS) has seen important developments [1- 
8], both in fundamental aspects and prospective tech¬ 
nological applications. As for the fundamental aspects, 
it was possible to understand the mechanisms of inter¬ 
action between diluted magnetic impurities allowing fer¬ 
romagnetic semiconductors at room temperature [7, 8]. 
From the point of view of prospective applications, the 
marriage between the world of two-dimensional semicon¬ 
ductors and the world of magnetic data storage could 
provide us two-dimensional spintronics devices such as 
non-volatile magnetoresistive memories, spin valve, spin- 
based transistors and even magnetically enhanced opto¬ 
electronics devices [9]. 

Theoretical works based upon first-principles calcula¬ 
tions have greatly contributed to the understanding of 
magnetism in dilute semiconductor materials [4, 5]. Fur¬ 
thermore, it has been possible to achieve a spontaneous 
magnetic order above room temperature in a few classes 
of materials, including magnetic semiconducting oxides 
and nitrides [G, 8]. However, in practice disorder effects 
and thermal fluctuations result in Curie temperatures 
well below the expectations. 

In parallel, despite the success of two-dimensional (2D) 
materials such as graphene [10, 11], transition metal 
dichalcogenides (TMD) [12, 13] and phosphorene, mag¬ 
netism in 2D semiconductors has remained almost unex¬ 
plored. One of the reasons for this might have been a 
theoretical result [14] indicating that the Curie temper¬ 
ature in 2D DMS tends to be substantially lower than 
for the corresponding 3D system. However, a subse¬ 
quent work dissipated doubts by showing that including 
an adequate description of the disorder and the temper¬ 
ature dependence of the Fermi energy, it is possible to 
achieve much higher transition temperatures [15]. Con- 
cordantly, a Monte Carlo study of Mn-doped M 0 S 2 , with 
first-principles parameters, has unraveled the potential 


of monolayer materials by predicting Curie temperatures 
above room temperature [16]. 

From a technical point of view, 2D semiconductors 
have other winning factors that can be explored in mag¬ 
netic or spintronic devices. Firstly, the carrier concen¬ 
tration can be controlled externally by gating. And 
there is room for improving the control over the impu¬ 
rity concentration, for example exploring the possibility 
of using adatoms to incorporate impurities in concen¬ 
trations above the solubility limit. In practice, stud¬ 
ies in magnetic semiconductor nanostructures with lower 
dimensionalities, including semiconductor nanocrystals 
and nanowires [6, 17-19] doped with transition metals 
(TM), indicate that the confinement effect and improved 
control of magnetic dopants can be used to increase the 
Curie temperature. 

Phosphorene [20-22], a monolayer of black phospho¬ 
rus, presents some advantages compared with other pre¬ 
viously studied 2D semiconductors. Different from semi¬ 
conducting TMD, it is an elemental semiconductor and 
therefore less prone to vacancies [23]. Further, it has an 
anisotropic structure allowing in principle for different 
exchange interaction along the two principal directions. 
In host heteropolar semiconductors as GaAs (zincblende 
structure), there is antisite defects which reduces the 
Curie temperatures of TM-doped GaAs dilute magnetic 
semiconductors [4]. It is expected that in 2D heteropo¬ 
lar semiconductors as TMD, antisite defects [24] will also 
cause this effect. However, in elemental host semiconduc¬ 
tors as phosphorene, this kind of defects does not exist. 
Thus, Curie temperatures of dilute magnetic semiconduc¬ 
tors based on phosphorene are likely to be higher than 
those found in TM-doped TMD. 

In this article, we use Co as a prototype magnetic im¬ 
purity to show that it possible to achieve magnetism 
(ferromagnetism and antiferromagnetism) in phospho¬ 
rene and that this can be switchable by the gate voltage 
through the control of the carrier density. We first stud¬ 
ied the stability and electronic properties of Co defects in 
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FIG. 1. Ball-and-stick representation of phosphorene with: (a) substitutional Co (Cop) defect and (b) interstitial (or adsorbed) 
Co (Coi) defect. Electronic band structures for: (c) pristine phosphorene, (d) Cop-doped phosphorene, (e)-(f) Coi-doped 
phosphorene for spin up and down, respectively. 


phosphorene (interstitial and substitutional), and there¬ 
after, the magnetic orders acquired through the long- 
range interaction between diluted Co defects. 


METHODS 

Electronic, energetic and magnetic properties of the 
systems studied were calculated via first principles meth¬ 
ods based upon density functional theory (DFT) [25, 
26]. The simulations were performed with the Vasp 
code (Vienna Ab initio Simulation Package) [27, 28] 
with external potentials in PAW approximation [29] 
and exchange-correlation functional in Perdew-Burke- 
Ernzerhof approximation (GGA-PBE) [30] and Heyd- 
Scuseria-Ernzerhof approximation (HSE06) [31]. We 
used cutoff energy of 400 eV for Kohn-Sham orbitals 
and first Brillouin zone /c-points sampling with 5x5x1 
grid calculated with Monkhorst-Pack algorithm [32] for 
5x4 supercell. For Density of States (DoS) and Pro¬ 
jected Density of States (PDoS), we use Brillouin zone 
integrations with fc-points sampling in 15 x 15 x 1 grid. 
The geometries were relaxed with residual force criteria 


of 10“^ eV/A. The primitive unit cell for pristine phos¬ 
phorene in our calculations was characterized by lattice 
constants a « 3.298 A and b ft! 4.625 A. We used vac¬ 
uum distances of 15 A to avoid spurious interactions. For 
the calculation of magnetic anisotropy energies (MAE), 
we used noncollinear spin polarization in fully-relativistic 
framework. Volumetric data were analyzed with Vesta 
program [33]. 

Ionisation levels were calculated using the marker 
method [34], taking the pristine supercell as reference. 
This method was found to give similar results to ionisa¬ 
tion levels calculated using the formation energy method, 
after correcting for spurious electrostatic interactions, 
which partially cancel when systems with the same 
charge are compared [35]. 


COBALT POINT DEFECTS 

Cobalt atoms can either be adsorbed on phosphorene 
or occupy a substitutional position (Fig. 1). Both of 
them lead to very small lattice distortion. Substitu¬ 
tional cobalt (Cop) forms three bonds, 3-5% shorter than 
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FIG. 2. (color online) (a)-(e) Density of States (DoS) shown as light blue and light red areas, and Projected Density of States 
(PDoS) on cobalt atoms shown as dark blue and dark red areas. The PDoS are decomposed on atomic orbitals: (a) ci^ 2 , (b) 
dxy, (c) d^ 2 _y 2 , (d) dxz, and (e) dyz. The DoS scale is ten times larger than the PDoS scale, (f) Schematic representation of 
3d level splitting in crystal field with Cs point symmetry and exchange interaction. 


the original P-P bonds [Fig. 1(a)]. Interstitial (or ad¬ 
sorbed) cobalt (Coi) also forms three bonds, connecting 
P atoms in two zigzag ridges on the same side of the layer 
[Fig. 1(b)]. The adsorption energy of Coi defect is among 
the highest of TM adatoms on phosphorene [36] . 

In order to calculate the electronic properties of an 
isolated Co point defect in phosphorene, were used a 
5x4 supercell, so that these Co point defects are di¬ 
luted in concentrations of 1.23% (interstitial) and 1.25% 
(substitutional). Unless otherwise stated, all the calcula¬ 
tions described in this article are in these concentrations. 
Whereas the Cop-doped phosphorene does not presents 
nonzero magnetic moments, the Coi-doped phosphorene 
presents a total magnetic moment of m = 1 /ip- This 
is 0.12 eV lower in energy than the nonmagnetic state 
(NM). 

The Magnetic Anisotropy Energy (MAE) for the di¬ 
luted Coi defect was obtained from fully-relativistic and 
noncollinear spin-polarised calculations. The vertical 
spin-polarization (z-direction) proved to be the lowest 
energy, whereas spin polarizations in other directions {x 
and y-direction) are related by MAE(?/z) = 1.0 meV, 
MAE(a; 2 ) = 0.8 meV, and MAE(a;j/) = 0.2 meV. Since 
the magnetic anisotropy is small, we will treat the mag¬ 
netic impurities using a collinear spin-polarised density- 
functional formalism. 


Electronic band structures 

The band structures calculated for pristine, Cop-doped 
and Coi-doped phosphorene are shown in Fig. l(c)-(f). 
For the spin-degenerated Cop-doped phosphorene band 
structure [shown in Fig. 1(d)], defect levels with disper¬ 
sion of approximately 0.2 eV were introduced at about 
mid-gap, and the valence band shifts in such way that 
the valence band maximum (VBM) becomes located at 
A-point. This shifts in VBM will suppress the lumines¬ 
cence of phosphorene [21, 37]. As for the band structures 
(spin up and down) of Coi-doped phosphorene, two local¬ 
ized levels are introduced in the spin down channel. One 
is occupied (below the Fermi level) and other unoccupied 
(above the Fermi level). For the spin up channel, all Co 
d-states are occupied, leading to an imbalance resulting 
in a local magnetic moment of 1 /xp, which gives rise to 
a macroscopic magnetisation. 

Origin of the defect state 

The origin of the local magnetic moment in Coi-doped 
phosphorene can be understood by analyzing the crystal 
field splitting, the hybridisation between the Co valence 
states (4s^3d^) and the phosphorene host states, which 
have mainly s and p character. 

From the symmetry breaking point of view, while pris¬ 
tine phosphorene has C 2 v point symmetry, cobalt point 
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defects reduces the symmetry to Cs ■ Only a mirror plane, 
passing through the defect and perpendicular to phos- 
phorene is preserved. In this lower symmetry, cobalt d 
orbitals transform as 3^' 0 2A”, where d^y, dx 2 _y 2 and 
dz 2 belong to the A' representation, and dxz and dyz be¬ 
long to the A" representation [38] . Phosphorus p orbitals 
splits in 2A' 0 A", where Px and Py are in A! represen¬ 
tation, and Pz is in A!' representation. We observe that 
the pd hybridisation occurs in A with p^, dxy and 
dx 2 -y 2 orbitals, and in A" with pz and dxz- 

Different from interstitial transition metals in other 
materials, the hybridisation is quite strong and therefore 
the electronic states cannot be understood based only 
on a crystal field splitting model. This is apparent in 
the Projected Density of States (PDoS) on cobalt atoms, 
obtained from the decomposition on atom-centred real 
spherical harmonics, shown in Fig. 2(a)-(e). We noticed 
that there is pd hybridisation between p orbitals from 
phosphorene and d orbitals from cobalt defects, resulting 
in bonding and antibonding hybrid states. The exchange 
field splits each of these states, leading to a splitting 
into four peaks for each d orbital. This effect is more 
pronounced on dxy, dxz and dx 2 -y 2 orbitals, which have 
maximum amplitude in the xy plane [Fig. 2(b),(c) and 
(d)]. The dz 2 and dyz orbitals are nodal at the xy plane 
and therefore undergo smaller hybridisation and have 
lower energies. A schematic representation of crystal field 
and exchange splitting is shown in Fig. 2(f). For simplic¬ 
ity, this takes into account only the states with highest 
d character for each d orbital type and spin. Whereas 
the dxz, dyz, dxy and dj .2 orbitals present exchange split¬ 
ting in the range 0.32-0.45 eV, the dx 2 -y 2 orbital present 
exchange splitting of approximately 1.08 eV. This large 
apparent exchange splitting is in fact a combined effect 
of pd hybridisation and exchange field splitting, resulting 
in higher peak for spin up bonding state and spin down 
antibonding state. 

Formation energy 

We define the formation energy for neutral defects 
(AO) as [39] 

Ef[X°] = Atotp :A°] - Etot[P] - Px - nppp, (1) 

where AtotP :A°] is the total energy of phosphorene with 
A° defect, i?tot[P] is the total energy of pristine phos¬ 
phorene, Px = PCo is the chemical potential of cobalt 
calculated from 3D crystal with hexagonal lattice, /ip 
is the chemical potential of phosphorus calculated from 
pristine phosphorene. The number of phosphorus atoms 
missing (np) is np = 0 for Coi defects and np = 1 for 
Cop defects. 

In the neutral charge state, substitutional cobalt de¬ 
fects have the lowest formation energy. While formation 
energy of neutral Cop is A'^[Cop] = 1.20 eV, neutral 



FIG. 3. (color online) Formation energy vs Fermi level for 
charged and neutral states of Cop (light green) and Coi (dark 
violet) defects. 

Coi have i?-^[Co?] = 1.42 eV. However, both defects are 
electrically active and therefore the formation energy de¬ 
pends on the Fermi level Ep (Fig. 3). Calculating the 
thermodynamic ionisation levels /(t/O) and 7(0/—) with 
the marker method [34], we can evaluate the formation 
energies of charged point defects as a function of the 
Fermi level as 

Ef[X+]iEF) = Ef[X°] - 7(0/0) 0 Ep, (2) 

for positive charged point defects (A+), and 

E^[X-]{Ep) = Af0 E^ [A°] - 7(0/-) - Ep (3) 

for negative charged point defects (A“), where the donor 
and acceptor levels are given relative to the valence and 
conduction band edges, respectively. The Fermi level 
is set to zero at VBM. In equation (3), is the 

bandgap at PBE level. Substitutional Co is an acceptor 
with 7(0/—) level at midgap (at the PBE level) and a 
double acceptor at Ep = 0.78 eV, whereas Coi is ampho¬ 
teric. Thus, in p-type material Coi is stabilised, becom¬ 
ing more energetically favourable than Cop. According 
to the PBE calculations, this happens for Ep between 7?„ 
and Ey + 0A6 eV. However, since the bandgap is severely 
underestimated at the PBE level, we extrapolated the 
results using the valence and conduction band shifts ob¬ 
tained from HSE calculations for the pristine phospho¬ 
rene. Here, we assume that the gap levels of Coi and 
Cop are very localised, are not shifted with respect to the 
vacuum level by the exchange component [40]. By fixing 
the PBE and HSE vacuum level at same energy, we ob¬ 
serve the shifts AEy = E^^^ — E^^^ = —0.48 eV in the 
valence band, and AEc = E^^^ — E^^^ = 00.31 eV 
in the conduction band. The bandgap increase from 
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FIG. 4. (color online) Spin density isosurfaces of ±2-10“^ bolir“® for: (a) antiferromagnetic (AFM) order, and (b) ferromagnetic 
(FM) order. Positive values of isosurfaces are shown in blue, while negative values are shown in red. (c) Energy difference 
between FM and AFM orders as function of carrier densities (electrons and holes) for n-doped Co, (blue), p-doped Coj (red) 
n-doped Cop (green) defects. 


= 0.91 eV at PBE level to = 1.70 eV at 

HSE level. Because of this large shift in the valence band 
at HSE level, the formation energy of the donor Coi*" can 
be up to 0.42 eV smaller than the formation energy of 
CoO. 

MAGNETIC ORDER 

The pd hybridisation in Coi-doped phosphorene is 
characterized by long-range and anisotropic interactions 
between cobalt defects through the valence band, known 
as Ruderman-Kittel-Kasuya-Yosida (RKKY) interac¬ 
tions [41-43]. These interactions have oscillating ex¬ 
change coupling with the distance, leading to ferromag¬ 
netic and antiferromagnetic orders. 

In order to determine the magnetic order of the fun¬ 
damental state, we use a 4 x 4 supercell with four Coi 
defects in equivalent positions. In these systems, the Coi 
defects are diluted in a concentration of 5.88 %. In neu¬ 
tral system, the ferromagnetic (FM) order was more ener¬ 
getically favourable than antiferromagnetic (AFM) order 
by AEfm-afm = —24 meV. These magnetic orders are 
shown in Fig. 4(a)-(b) by spin density (s = — p^) iso¬ 


surfaces. Positive values of isosurfaces (s = -1-2 • 10“^ 
bohr“^) are shown in blue, and negative values (s = 
—2 • 10“^ bohr“^) are shown in red. 

Since Coi defects are amphoteric, they can accept or 
donate electrons. Codoping the system with electrons 
or holes, we can change the interaction between defects, 
and therefore, change the magnetic order of the sys¬ 
tem. We observe that when Coi-doped phosphorene is 
n-type, AEfm-afm decreases until a minimum value 
at 0.125e/Co (n = 1.7 • 10^^/cm^), and then increases 
to the ferromagnetic order is no longer stable at le/Co. 
For the p-type Coi-doped phosphorene, the difference 
AEpM-AFM increases monotonically with the density of 
holes, so that there is a phase transition from ferromag¬ 
netic to antiferromagnetic at a critical density of holes 
{pc = 1.11 • 10^^/cm^). 

For similar system, but with four Cop defects in phos¬ 
phorene, due to acceptor character of Cop defects, the 
magnetic orders were studied only codoping the system 
with electrons. In these systems, the Cop-doped phos¬ 
phorene begins from nonmagnetic order at neutral state 
to a ferromagnetic order with AEfm-afm = —19 meV 
at le/Co {n = 1.64-10^"^/cm^). Similar electric-field con¬ 
trol of magnetism was observed in (In,Mn)As quantum 
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well [2], 


CONCLUSION 

In summary, we show that interstitial and substitu¬ 
tional cobalt defects in phosphorene are stable, with low¬ 
est formation energy for Coi defect in p-doped phos¬ 
phorene, and Cop defect in neutral and n-doped phos¬ 
phorene. Also we show that through electron and hole 
codoping which can be controlled by gate potential, the 
exchange interaction between impurities can be engi¬ 
neered to yield ferromagnetic and antiferromagnetic or¬ 
ders. This effect can be used in prospective spintronic ap¬ 
plications. Whereas TM-doped M 0 S 2 have shown diluted 
magnetism with very high Curie temperatures (above 
room temperature) [16], dilute magnetic phosphorene 
presents some advantages due to absence of antisite de¬ 
fects and less prone to vacancies, indicating potential 
higher Curie temperatures and control of disordered im¬ 
purities. 
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